M^D candidate nucleus 05 Rh in relativistic mean-field approach 
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Abstract 

Following the reports of two pairs of chiral doublet bands observed in 105 Rh, the adiabatic and 
configuration-fixed constrained triaxial relativistic mean-field (RMF) calculations are performed 
to investigate their triaxial deformations with the corresponding configuration and the possible 
multiple chiral doublet (M%D) phenomenon. The existence of M%D phenomenon in 105 Rh is 
highly expected. 
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Chirality is a topic of general interest in nature science, such as chemistry, biology and 
physics. The occurrence of chirality in atomic nuclear structure was suggested for triaxi- 
ally deformed nuclei in 1997 [l| and the predicted patterns of spectra exhibiting chirality 
- chiral doublet bands — were experimentally observed in 2001 |2j. In addition to the 
triaxial deformation, the configuration with high-j valence particle(s) and valence hole(s) is 
also essential for chirality in nuclei. Up to now, candidate chiral doublet bands have been 
proposed in a number of odd-odd, odd-A or even-even nuclei in the A ~ 100, 130, 190 mass 
regions, for a review, see e.g. {3} Q|. 

Theory wise, chiral doublet bands were first investigated in the one-particle-one-hole- 



rotor model (PRM) and the corresponding tilted axis cranking (TAC) approximation 
Later on, numerous efforts have been devoted to the development of TAC fl-EJ and PRM 
models |8l-ll2| to describe chiral rotation in atomic nuclei. 

RMF theory has received wide attention due to its massive success in describing nuclear 



global properties and exotic phenomena I13H15I1. In order to describe the nuclear rotation 

nri 

phenomena, the cranked RMF theory in the context of principal axis rotation |16l . I17| as well 
as three-dimensional rotation 18] were developed. However due to numerical complexity, it 
was restricted to two-dimensional studies only, i.e., the magnetic rotation jisl. 

Based on the adiabatic and configuration-fixed constrained triaxial RMF calculation, 
triaxial shape coexistence with high-j proton hole- and neutron particle-configurations, are 
found in 106 Rh, which suggests the existence of a new phenomenon — M^D 20]. This 
prediction holds true for other rhodium isotopes 104 ' 106 ' 108 ' 110 Rh as well 21] . In particular, 
the prediction of the M%D in 106 Rh remains even with the time-odd fields indued 22]. 

Recently, two pairs of chiral doublet bands have been respectively observed in 105 Rh with 
three-quasiparticle configurations 71-519/2 (8> ^^11/2(^7/2,(^5/2) [23] and 7rg 9 / 2 <8> vh^^ [24]. It 
is interesting to verify whether these observations are the predicted MyD or not in the 



22] , the time-odd 



configuration-fixed constrained triaxial RMF approach. Similar as in Ref. 
fields are included. 

The starting point of the RMF theory is the standard effective Lagrangian density con- 
structed with the degrees of freedom associated with nucleon field (ip), two isoscalar me- 
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15] . Under 



son fields (a and u^), isovector meson field and photon field (A^) 
"mean-field" and "no-sea" approximations, one can derive the corresponding energy density 
functional, from which one finds immediately the equation of motion for a single-nucleon 



2 



orbital ipi(r) by variational principle, 

{a-\p- V(r)} + /3M*(r) + V (r)}^{r) = e^(r), (1) 

where M*(r) = M + g a a(r), M the mass of bare nucleon, Vo(r) the time-like component 
of vector potential, and V(r) the space-like components of vector fields. The details about 
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26). 



the solution of Dirac equation (|T|) with time-odd fields can be found in Refs. 

In the configuration-fixed constrained calculation, the same configuration is guaranteed 
during the procedure of constraint calculation with the help of "parallel-transport" J27!]. In 
addition to the /3 2 -constrained calculation 20| . the constraints on the axial and triaxial mass 
quadrupole moments are also performed to obtained the potential energy surfaces (PES) in 
the two-dimensional /S-7 plane 



In the present calculations, each Dirac spinor ipi(r) is expanded in terms of a set of 
three-dimensional harmonic oscillator (HO) basis in Cartesian coordinates with 12 major 
shells and the meson fields with 10 major shells. The pairing correlations are quenched by 
the unpaired valence nucleons in 105 Rh and thus neglected here. The effective interaction 
parameter set PK1 
into account by 



29j is applied. The center-of-mass (cm.) correction 29|, |30| is taken 

1 



£c m m.'- 2M ( P c.m.>> (2) 

where P c . m . is the total momentum operator of a nucleus with mass number A. In or- 
der to save time, the constrained triaxial RMF calculations without the time-odd fields 
are performed to search for the triaxial deformation parameters and the valence nucleons 
configuration favorable for chirality, which will be later confirmed by the calculations with 
time-odd fields. 

The potential energy surface in the f3-j plane (0^7^ 60°) for 105 Rh in the adia- 
batic constrained triaxial RMF calculations with PK1 is shown in Fig. [TJ All energies are 
normalized with respect to the binding energy of the absolute minimum, and the contours 
join points on the surface with the same energy (in MeV). The energy separation between 
contour lines is 0.4 MeV. From Fig. [U it shows that the ground state of 105 Rh is triaxially 
deformed with (5 = 0.28 and 7 = 23°. The second minimum is located at the area with 
(5 pd 0.4 and 7 ~ 5°. The energy difference between these two minima is about 0.5 MeV and 
corresponding barrier height is about 1.5 MeV. The behavior of shape coexistence is clearly 
shown. 



3 



PK1 




FIG. 1: (Color online) Contour plots of potential energy surface in /3-7 plane (0^7^ 60°) for 



29(. All energies 



105 Rh in constrained triaxial RMF calculations with effective interactions PK1 
are normalized with respect to the binding energy of the absolute minimum (in MeV). The energy 
separation between contour lines is 0.4 MeV. 
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FIG. 2: (Color online) The energy surfaces in adiabatic (open circles) and configuration-fixed (solid 
lines) constrained triaxial RMF calculation using effective interaction PK1 for 105 Rh. The minima 
in the energy surfaces for fixed configuration are represented as stars and labeled respectively as 
A, B, C, D, E, and F. Their corresponding triaxial deformation parameters f3 and 7 are also given. 

The /3 2 -constrained calculations, in which the triaxial deformation 7 is automatically 
obtained by minimizing the energy, have also been performed. In Fig. [21 the potential energy 
surfaces as a function of /3 in adiabatic (open circles) constrained triaxial RMF calculation 
with PK1 for 105 Rh are given, where four minima observed in the potential energy surfaces 
are labeled with A, B, C, and D respectively. 
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TABLE I: The total energies E tot , triaxial deformation parameters /3,7 as well as their correspond- 
ing valence nucleon configurations of minima for A-F in the configuration-fixed constrained triaxial 
RMF calculations. 



State 


Configuration 






Etot 




7 


Valence nucleons 


Unpaired nucleons 


A 


tt1c/- 3 2 <g> v±h\ l/2 






-898.09 


0.28 


24.7° 


B 


irlg-J 2 ®v(lg-J 2 2d\ /2 ) 




Trl^/2 


-897.76 


0.23 


21.3° 


C 


Tr(lg-p.p-j 2 lg^ /2 ) ® u{lg-^3s 2 1/2 lh 2 n/2 ) 




7r2p-^ 


-897.45 


0.41 


4.1° 


D 


7r ( 1 % 4 2 1 »7/2) ® K^T/Ml^) 




Trlg) /2 


-897.12 


0.37 


8.3° 


E 




vrl3 9 / 2 < 


% v{lh\ l/2 2d\ /2 ) 


-897.40 


0.25 


23.1° 


F 


7rl^ ® u(lh\ 1/2 lh\ 1/2 ) 


^9/2 5 


$ is(lh\ l/2 lh\ 1/2 ) 


-896.40 


0.28 


21.8° 



The total energies E tot , triaxial deformation parameters (3, 7 as well as their corresponding 
configurations of minima for A-D in the constrained triaxial RMF calculations are presented 
in Table [H Here, the state A represents the ground state, with triaxial deformation f3 = 0.28, 
7 = 24.7° and the corresponding valence nucleon configuration nlg~J 2 Cg> ulh^^. Note that 
the two of three proton holes in 1^9/2 orbital and two neutrons in the I/H1/2 orbital are 
pairwise, i.e., they occupy the degenerate time-reversal conjugate levels and don't contribute 
to the angular momentum. The corresponding unpaired nucleon configuration is 7ilg^ 2 . For 
state B, the triaxial deformation parameter is {(5 = 0.23, 7 = 21.3°), with corresponding 
valence nucleon configuration 7tlg~^ 2 ® u(lg~ 2 2 2d^ 2 ) and unpaired nucleon configuration 
7vlg^ 2 . Although A and B have different valence nucleon configurations, they have the same 
unpaired nucleon configuration. 

The triaxial deformation parameters for excited state C and D are (/? = 0.41, 7 = 
8.3°) and (/3 = 0.37, 7 = 8.3°), respectively. The valence nucleon configuration for C is 

fK^a^a/UflW ® K 1 07/2 3s i/2 1 ' l ii/2) and for D ^(^g/t^y ® v { l 9y 2 lh ii/2)i while the 
unpaired nucleon configuration n2p^j 2 for C and nlgy 2 for D. The quadrupole deformations 

of states C and D are larger than states A and B, which is due to the deformation driving 

orbital ^1^7/2 in states C and D. It should be pointed out that states A, B, C and D don't 

have the suitable particle-hole configurations for chirality. 



In the Refs. 



23 



the observed partners bands with the configurations irlg g ^ 2 
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vlh\ l / 2 (lg7/2, 2d 5 / 2 ) 1 and irlggh ® ^(1^11/2^11/2) are respectively suggested as candidate 
chiral doublet bands. It is interesting to examine the triaxial deformation of these states. 
For this purpose, the configuration-fixed constrained calculation with the suggested config- 
urations are performed and their corresponding energy surfaces are given in Fig. [2] in blue 
solid lines, and the energies, configurations and triaxial deformation parameters are listed 
in Table [H The excitation energy for minima E is 0.69 MeV and for F 1.69 MeV. Fur- 
thermore, the triaxial deformation suitable for chirality is found for E and F, namely 23.1° 
and 21.8° respectively, which together with the corresponding high- 7 proton hole and high-j 
neutron particle configurations will lead to the M%D phenomenon 20|. It should be noted 
that the present RMF calculations are restricted to the non-rotating mean field only, while 
the TAC calculations with Woods-Saxon potentials have demonstrated that the rotat ing 
mean field will become chiral for configuration E at fko = 0.20 MeV (strong paring) 23] 
and for configuration F at tko = 0.60 MeV (zero pairing) or 0.45 MeV (strong paring) {24]. 
Recently, the doublet bands with configuration F in 105 Rh are also investigated by triax- 
ial PRM calculations and the evolution of the chiral geometry with angular momentum is 
discussed [31]. 



The neutron and proton single-particle levels as functions of deformation j3 are given 
in Fig. [31 obtained in the adiabatic constrained triaxial RMF calculations for the region 
0.247 < (3 < 0.315. The occupations of neutron and proton for A, E and F are given in 
Fig. [3]^a) and Fig. [3]^b) respectively. The single-particle levels with positive (negative) parity 
are marked by solid (dashed) lines, and occupations are represented by filled circles (two 
particles) and stars (one particle). It should be pointed out that every single-particle level 
is degenerated with opposite simplex quantum number [l9|. It is easy to see that for ground 
state A, one valence proton occupies the 1*79/2 orbit and the last two neutrons occupy the 
degenerated time reversal conjugate orbit l/in/2- For the positive parity state F, the two 
I/111/2 neutrons are unpaired (in contrast to the state A), i.e., they occupy two different 
l/in/2 levels. For the negative parity state E, the two unpaired neutrons occupy the l/iu/2 
and 2g?5/2 levels. It is interesting to note that states E and F compete with each other in 
energy. However, due to different parities, states E and F do not mix up and could contribute 



the M%D phenomenon 



20]. 



Similar as in Ref . 22[ , the calculations with time-odd fields are also performed to confirm 
the above discussions. In Table HT[ the calculated triaxial deformation parameters (3, 7, total 
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FIG. 3: (Color online) Neutron and proton single-particle levels obtained in constrained triaxial 
RMF calculations with PK1 as functions of deformation f3. Positive (negative) parity states are 
marked by solid (dashed) lines. Occupations corresponding to the minima in Fig. [2] are represented 
by filled circles (two particles) and stars (one particle). 
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TABLE II: The triaxial deformation parameters /?, 7, total energies E tot and the excitation energies 
E x for states E and F calculated with and without time-odd fields, compared with the experimental 
bandhead energies of two rotational bands based on the configurations of E and F. The spin and 
parity P r for ground state and band head of two rotational bands are also given. 



State 


Time- 


■even 


Time-odd 


exp 


(A 7) 


-E-x(-Etot) 


(0.7) 


E x {E to t) 


r e x 


A 


(0.28,24.7°) 


(-898.09) 


(0.28,24.9°) 


(-898.14) 


7+ 
2 


E 


(0.25,23.1°) 


0.69(-897.40) 


(0.25,23.3°) 


0.37(-897.76) 


^~ 2.417 


F 


(0.28,21.8°) 


1.69(-896.40) 


(0.28,22.1°) 


1.43(-896.70) 


f + 2.982 



energies E tot , and the excitation energies E x for states E and F with and without time-odd 
fields are given and compared with the experimental bandhead energies of two rotational 
bands. The effects of the time-odd fields on the triaxial deformation parameters /3,7 are 
negligible. Their contribution for the total energy are respectively —0.05 MeV for state A, 
-0.3 MeV for state E and -0.36 MeV for state F. 

The experimental spin, parity and excitation energies of bandheads for two rotational 
bands based on the configurations of E and F are , 2.417 MeV and f + , 2.982 MeV 
respectively. It can be seen that the displacement for the energies of two bandheads (0.57 
MeV) has been reasonable reproduced by the RMF calculations (1.0 and 1.06 MeV). 
It is necessary to note here that in 105 Rh, the suggested configuration 7rlg~i_ 



vlh\ 1 / 2 (lg 7 /2,2d 5 /2) 1 for the candidate chiral doublet bands (band 7 and 8 in Ref. [231 ]) 



involves the orbits of pseudospin doublet states (1^7/2, 2^5/9). So a competing interpretation 
of band 7 and 8 includes the pseudospin doublet bands [3j. Further efforts are needed to 
address this point. 

In summary, adiabatic and configuration-fixed constrained triaxial RMF approaches have 
been applied to investigate the M%D candidate nucleus 105 Rh. According to the suggested 
high-j proton hole and high-j neutron particle configurations, their triaxial deformations fa- 
vorable for the construction of the chiral doublet bands are obtained from the configuration- 
fixed constrained triaxial RMF calculations. The existence of M%D phenomenon is expected 
in 105 Rh. Future efforts should be made to investigate whether the rotating mean field will 
attain chirality or not. 
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